Structure of Complex-Periodic and Chaotic Media with Spiral Waves 
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The spatiotemporal structure of reactive media supporting a solitary spiral wave is studied for 
systems where the local reaction law exhibits a period-doubling cascade to chaos. This structure is 
considerably more complex than that of simple period- 1 oscillatory media. As one moves from the 
core of the spiral wave the local dynamics takes the form of perturbed, period-doubled orbits whose 
character varies with spatial location relative to the core. An important feature of these media is 
the existence of a curve where the local dynamics is effectively period- 1. This curve arises as a 
consequence of the necessity to reconcile the conflict between the global topological organization 
of the medium induced by the presence of a spiral wave and the topological phase space structure 
of local orbits determined by the reaction rate law. Due to their general topological nature, the 
phenomena described here should be observable in a broad class of systems with complex periodic 
behavior. 
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I. INTRODUCTION 

Spiral waves are spatio-temporal patterns typically 
found in distributed media with active elements. They 
have been studied extensively for excitable and oscilla- 
tory media. For both types of media, it is conven- 
tional to consider systems with two dynamical variables. 
Activator-inhibitor or propagator-controller systems are 
often used to analyse spiral dynamics in excitable me- 
dia while the complex Ginzburg-Landau equation 
is the prototypical model describing spatially-distributed 
oscillatory media near the Hopf bifurcation point Q] . 

Spiral waves may also exist in media where the lo- 
cal dynamics supports complex periodic or even chaotic 
motion that cannot be represented in a two-dimensional 
phase plane. Various patterns involving rotating spiral 
waves have been observed in coupled map lattices or 
reaction-diffusion dynamics based on the Rossler chaotic 
attractor ||. The three- variable reaction-diffusion sys- 
tem with chaotic local reaction kinetics given by the 
Willamowski- Rossler rate law j(| has been studied in . 
Stable spiral waves exist in this system and the nucleation 
and annihilation of spiral pairs leading to spiral turbu- 
lence have been observed. The change of dimensionality 
of phase space from two to three significantly complicates 
the description of the dynamics. Descriptions in terms 
of phase and amplitude, well established for two- variable 
models, cannot be directly generalized. Although several 
definitions have been proposed for the phase of chaotic 
oscillations, all of them suffer from some degree of ambi- 
guity (see H for a discussion). Similar difficulties arise 
in the consideration of nonchaotic oscillatory dynamics 
which is nevertheless more complex than a single loop in 
phase space; for example, in the oscillations that appear 
in the period doubling cascade to chaos or in the mixed- 
mode oscillations observed in experiments in chemical 
systems. H 

In this paper we study the spatiotemporal organiza- 



tion of a reacting medium which supports a single spiral 
wave and where the local rate law exhibits period-2™ or 
chaotic oscillations. Through an analysis of the dynamics 
at different spatial points in the medium we show that a 
number of phenomena arise for n > which are nonexis- 
tent in period- 1 oscillatory media. Section [n] introduces 
the model and presents some features of the spiral wave 
behavior in a chaotic medium. The local dynamics in the 
medium is considered in detail in Sec. III. The analysis 



allows one to identify the loop exchange process for local 
trajectories and the complicated pattern of the distribu- 
tion of different types of local dynamics in the medium. A 
characteristic feature of this distribution is the existence 
of a curve where the local dynamics is effectively period- 
1. Section IV introduces a coarse-grained description of 



2™-periodic local orbits which allows one to character- 
ize the local dynamics that is observed in the medium. 
The topological conflict between the phase space struc- 
ture of local trajectories and the constraints imposed on 
the medium by the existence of a spiral wave is consid- 
ered in Sec. [V| We show that the observed changes of 
the local orbits are necessary to maintain the global co- 
herence of the medium. The conclusions of the study are 
presented in Sec. VI. 



II. 



SPIRAL WAVES IN PERIODIC AND 
CHAOTIC MEDIA 



While many aspects of the phenomena we describe in 
this paper are general and apply to systems in which com- 
plex periodic or chaotic orbits exit, we consider situations 
where a chaotic attractor arises by a period-doubling cas- 
cade and confine our simulations to the Willamowski- 
Rossler (WR) model [§, 

Ai + X ^ 2A, X + Y ^ 2V, 
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Only the X, Y and Z species vary with time; all oth- 
ers are assumed fixed by flows of reagents. Study of this 
model allows us to illustrate most features of the struc- 
ture of a spatially distributed medium supporting spiral 
waves. In addition, it is useful to deal with a specific ex- 
ample since certain aspects of the analysis of periodic and 
chaotic orbits in high-dimensional concentration phase 
spaces rely on geometrical constructions that pertain to 
a specific class of attractors. 

The rate law that follows from the mechanism (fiT) is 



dc x {t) 
dt 

dc yW) 
dt 

dc z {t) 
dt 



KlC x - K~lC x - K 2 C x C y + K-%<? v - KiC x C z 
+ K_ 4 = R x (c(t)) , 



K 2 C x C y - K- 2 C y - K 3 C y + K-3 



(2) 



R y (c(t)) 



-K4C X C Z + K-4 + K 5 C Z - K_ 5 Cg = R y (c{t)) , 



where the rate coefficients k, include the concentrations 
of any species held fixed by constraints. We take k 2 
to be the bifurcation parameter while all other coeffi- 
cients are fixed: (k\ = 31.2, K—\ — 0.2, K- 2 = 0.1, K3 = 
10.8, K_3 = 0.12, k 4 = 1.02, k_ 4 = 0.01, k 5 = 16.5, k_ 5 = 
0.5). In this parameter region the WR model has been 
shown jlO) to possess a chaotic attractor arising from a 
period-doubling cascade as n 2 is varied in the interval 
[1.251,1.699]. 




FIG. 1. Chaotic attractor for the Willamowski-Rossler 
model at K2 = 1.567. 

Figure [j] shows the four-banded chaotic attractor at 
k 2 = 1.567. Throughout the entire parameter domain 
k 2 £ [1.251,1.699] the system's attractor is oriented so 
that its projection onto the (c x , c y ) plane exhibits a folded 
phase space flow circulating around the unstable focus c* . 
This allows one to introduce a coordinate system in the 
Cartesian (c x , c y ,c z ) phase space which is appropriate for 



the description of the attractor. We take the origin of a 
cylindrical coordinate system (p, 4>, z) at c* so that the z 
and zero-phase-angle (4> = 0) axes are directed along the 
c z and c v axes, respectively. The phase angle cj> increases 
along the direction of flow as shown in Fig. 0. 




FIG. 2. Cylindrical coordinate frame (p, cf>, z) with origin at 
c* in the (c x ,Cy,c z ) phase space. A period-2 orbit is shown 
in this coordinate frame. 

For a period- 1 oscillation <f> coincides with the usual 
definition of the phase and uniquely parametrizes the at- 
tractor p a — p a {4>),z a = z a {(j>),(j) £ [0, 2tt). After the 
first period-doubling this parametrization is no longer 
unique since the periodic orbit does not close on itself 
after </> changes by 27r. For a period-2™ orbit 2™ of its 
points lie in any semi-plane (f> = (f>o. The angle variable 
$ € [0, 2™ • 2ir) may be used to parametrize the period-2™ 
attractor if one acknowledges that all $ from the interval 
[0, 2" • 2ir) are different but any two values of $, $1 and 
$2, with $2 = $1 + 2™ • 27T, are equivalent. For a chaotic 
orbit (n — » 00) all angles $ € [0, 00) are non-degenerate. 
When $ is defined in this way it is no longer an observ- 
able. Indeed, any $ £ [0,2™ • 27r) can be represented as 
<f> = 4> + m ■ 2n where (f> £ [0, 2 7r) and m £ N. While 
4> is just the angle coordinate in the (p, 0, z) system and 
is a single-valued function of the instantaneous concen- 
trations <f) = cj){c x {t),Cy{t),c z {t)), the integer number of 
turns to can be calculated only if the entire attractor is 
known. 

The spatially-distributed system is described by the 
reaction-diffusion equation, 



dc(x,t) 
dt 



R(c(x,i)) +L>V 2 c(x,t) , 



(3) 



where we have assumed the diffusion coefficients of all 
species are equal. If the rate law parameters correspond 
to a period-1 limit cycle, we may initiate a spiral wave 
in the medium and describe its dynamics and structure 
using well-developed methods. The core of such a spiral 
wave is a topological defect which is characterized by the 
topological charge [jllj 



^ I V0(r) ■ dl = m 



(4) 
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where <p(r) is the local phase and the integral is taken 
along a closed curve surrounding the defect. To obtain 
additional insight into the organization of the medium 
around the defect the local dynamics may be considered. 
For this purpose we introduce a polar coordinate system 
r = x r<j(i) = (r, 6) centered at the defect whose (pos- 
sibly time-dependent position) is r<z(i). Let c(r, t) be a 
vector of local concentrations at space point r = (r, 8) . 
A local trajectory in the concentration phase space from 
t — to to t = to + t at point r in the medium will be 
denoted by 

C(r|to,T) = {c(r,t)|t G [to,t + T]} . (5) 

Figure [| shows a number of local trajectories C(r, 9\to, r) 
at points with increasing separation r from the defect for 
a period-1 oscillation at n 2 = 1.420. One sees that as 
r — » the oscillation amplitude decreases and the limit 
cycle shrinks to the phase space point c % corresponding 
to the spiral core. The results of our simulations show 
that the value of c* d differs only slightly from c* which 
is chosen as the origin of the coordinate frame (p,<f>,z). 
Thus, the angle 4> can serve as a phase that characterizes 
all points in the period-1 oscillatory medium except for a 
small neighborhood of the defect with radius r « 1. Jl^] 
The concentration field c(r, t) is organized so that the in- 
stantaneous (c x , c y ,c z ) phase space representation of the 
local concentration on any closed path in the medium 
surrounding the defect is a simple closed curve encircling 
c*. For large r, r > r rnax (in Fig. | r max rj 40), one finds 
that C(r,0\to,T) ceases to change shape and is indistin- 
guishable from the period-1 attractor of (||) on the scale 
of the figure. 




FIG. 3. Local trajectories calculated for the period-1 oscil- 
latory medium (« 2 = 1.420) at radii 5, 10, 20, 30, 40, 56, and 
fixed 8. The periodic orbits grow monotonically in size with 
r; the difference between trajectories corresponding to r — 40 
and r — 56 is not resolved on the scale of the figure. Local 
orbits appear to be independent on angle 9. The location of 
c* is designated by a diamond. 

One may initiate the analog of a defect in 2™-periodic 
and chaotic media. The defect serves as the core of a 
spiral wave which may exist even if the oscillation is not 
simply period-1. A defect was introduced in the center of 



the medium by fixing c z (r) = c* and choosing initial con- 
centrations (c x (r), Cy(r)) to produce orthogonal spatial 
gradients. The influence of the symmetry of the spatial 
domain on the dynamics was investigated by performing 
simulations on square (L x L) arrays as well as on disk- 
shaped domains with radius R. No-flux boundary condi- 
tions were used to prevent the formation of defects with 
opposite topological charge within the medium and to 
minimize effects arising from the self-interaction of spiral 
waves. The implementation of these initial and boundary 
conditions does not guarantee the formation of a solitary 
stable spiral wave; new spiral pairs and other patterns 
(e.g. pacemakers) may appear as a result of instabilities 
of the spiral arm and lead to spiral turbulence. The abil- 
ity to maintain a stationary spiral wave in the center of 
the medium is sensitive to the parameters. For various 
values of the system size and rate constants the defect 
can move along expanding or contracting spiral trajecto- 
ries or trajectories with complex "daisy" -like forms [jL3j. 
Simulations show that the stability of a spiral wave with 
a stationary core located at the center of the medium in- 
creases with the system size and for rate constants lying 
close to the chaotic regime within the period-doubling 
cascade. In the following we restrict our considerations 
to parameters that lead to the formation of a single spiral 
wave whose core is stationary and lies in the center of the 
domain. Long transient times (« 10 2 spiral revolutions) 
are often necessary to reach this attracting state. 




FIG. 4. Frames showing a rotating spiral wave in the 
chaotic (ft2 = 1.567) disk-shaped medium with R — 80. The 
local angle variable 4>{r, 6, t) is shown as grey shades. Time 
increases from left to right and from top to bottom. The 
frames are separated by one period of spiral revolution T r . 
The integration time step is At = 10 -4 and the scaled diffu- 
sion coefficient is DAt/(Ax) 2 = 10~ 2 . 

Figure |^ shows four consecutive states of the disk- 
shaped medium with R — 80, separated by one period 
of the spiral rotation, T r , for n 2 — 1.567 where the rate 
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law supports a chaotic attractor. Only within a suffi- 
ciently small region with radius r k 20 centered on the 
defect does the medium return to the same state after one 
period of spiral rotation. At points farther from the de- 
fect the system appears to return to the same state only 
after two spiral rotation periods. The transition from 
period- 1 to period-2 behavior occurs smoothly along any 
ray emanating from the defect. 

III. ANALYSIS OF LOCAL DYNAMICS 

More detailed information may be obtained from an 
investigation of the local dynamics of the medium sup- 
porting a spiral wave. Local trajectories C(r\ta,T) were 
computed along rays emanating from the defect at var- 
ious angles 8. Figure || (left column) shows short-time 
trajectories (t w 10T r ) at different radii r and arbitrary 
but large to. These trajectories clearly demonstrate that 
the local dynamics undergoes transformation from small- 
amplitude period- 1 oscillations in the neighborhood of 
the defect to period-4 oscillations near the boundary. |l~4| ] 




FIG. 5. Local trajectories C(r\to,r) for the disk-shaped 
medium (k 2 = 1.567, R = 80): (a, d) r = 10; (b, e) r = 35; 
(c, f) r = 76. The observation times are r « 10T r for left 
column and r = 80T r for right column. All the trajectories 
are shown on the same scale. 

The well-resolved period-doubling structure of 
C(r\to,r) is destroyed if the time of observation r be- 
comes sufficiently large. The right column of Fig. shows 
trajectories sampled at the same spatial locations but 
with the time of observation r = 80T r . These long-time 
trajectories appear to be "noisy" period-1 and period-2 
orbits: the trajectory in panel (d) is a thickened period-1 
orbit while both the period-2 (panel (b)) and period-4 
(panel (c)) orbits now appear as thickened period-2 or- 
bits in panels (e) and (f) with trajectory segments lying 
between the period-2 bands. As r tends to infinity the 
resulting local attractor, C(r), is independent of to and 
the angle 9. 



A. First-return maps 



An analysis of the local trajectories shows that the 
period-doubling phenomenon is not a monotonic func- 
tion of r. Consider the first return map constructed from 
a Poincare section of a local trajectory C(r\to,r) in the 
following way: choose the plane c y = c* with normal 
n along the c y axis as the surface of section and select 
those intersection points where n forms a positive angle 
with the flow. This yields a set {(c x (r,t n ),c z (r,t n ))\n € 
[1,7V]} where t < t x < t 2 < ... < t N < t + r 
is a sequence of times at which the trajectory crosses 
the surface of section. For the WR model the points 
(c x (r, t n ), c z (r, t n )) lie on a curve which deviates only 
slightly from a straight line. Consequently, we may 
choose either c x or c z to construct the first return map. 
Let £n(r) = c x (r, t n ) denote a point in the Poincare sec- 
tion. The relation £ n +i(r) = /(£ n (r)) between the suc- 
cessive intersections of the Poincare surface defines the 
local first return map, 

fl(r|t 0l r) = {(£n(rU„+i(r)) 

\t n £ [t ,t + T],n 6 [l,N]} . (6) 

Combining such maps for all r along some ray emanating 
from the defect at an angle 6, we obtain the cumulative 
first return map, 

G(e\t ,r)= |J g(r,6\t ,T). (7) 
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FIG. 6. Cumulative first return map G constructed for the 
disk-shaped array (ki — 1.567, R = 80). The letters indicate 
the r values discussed in the text: (a) 20, (b) 31 and (c) 43. 
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For sufficiently long times r, g is independent of 9 
and to. Letting linv^oo g(r, 9\to, r) = g(r), we may 
write the corresponding cumulative first return map as 
G = lim r _, 00 G(9\t ,r). Figure || shows G for the disk- 
shaped medium under consideration. The first return 
map is comprised of several branches which can be iden- 
tified as thread-like maxima of the first return map point 
density. These branches are parametrized by the spa- 
tial coordinate with r increasing from the bottom left 
corner to the ends of the wide-spread arms of G (cf. 
Fig. ||). Generally for r < 40 points lying on lines 
+ S,n+i(r) = const belong to the same g(r) though 
overlaps of neighboring <?-map points are common. Thus, 
measuring the separation between branches of G in the 
direction perpendicular to the bisectrix one can deter- 
mine the character of C(r). In spite of some evidence 
of fine structure, from the fact that map points are lo- 
cated along the bisectrix in Fig. || one can infer that up 
to r = 20 the local dynamics is predominantly period- 1. 
Starting from r = 21 (labeled by a in Fig. |), G splits 
into two branches which diverge from the bisectrix indi- 
cating a period-2 structure of C(r). As r increases these 
branches bend and cross the bisectrix at r = 31 (labeled 
by b in Fig. ^J), indicating a return of the local dynam- 
ics to the period- 1-like pattern. After this crossing the 
separation between the branches grows rapidly reflecting 
the development of period-2 structure. An examination 
of the main branches of G reveals period-4 fine structure. 
This period-4 structure is visible for r > 28 and beyond 
r « 43 (labeled by c in Fig. ^|) it becomes prominent and 
can be easily seen in the structure of C(r) (cf. Fig. |^). 



B. Loop exchange and Q curve 

From the analysis of the time series of the local concen- 
tration one may determine the processes responsible for 
the differences between the local trajectories C(r\to,r) 
for short and long time intervals r (cf. Fig. ||). Fig- 
ure shows the signature of this phenomenon for c x (r, t) 
at r = 50 in a disk-shaped array with R = 80 and 
«2 = 1.544, a parameter value corresponding to period- 
4 dynamics in the rate law. Every second maximum of 
c x (r,t) is indicated by diamond or cross symbols. The 
envelope curves obtained by joining like symbols cross at 
t — t ex , thus the curve which connected large-amplitude 
maxima at t < t ex joins low-amplitude maxima at t > t ex 
and vice-versa. This implies that if at some to < t ex 
the representative point c(r, to) was found on the small- 
amplitude band of period-2, then at t = to + nT2 > t ex , 
where Ti is the period of the period-2 oscillation, it will 
be found on the larger-amplitude band. p"5| This phe- 
nomenon can be interpreted as an exchange of the lo- 
cal attractor's bands. Indeed, approaching t ex from the 
left one finds that with each period of oscillation the 
small-amplitude band grows while large-amplitude band 
shrinks. At t = t ex both bands reach and pass each other. 




FIG. 7. Concentration time series c x (r,t) at r — 50 for the 
disk-shaped array (/t2 = 1.544, R = 80) showing the loop 
exchange process. Time unit equals 10 At. 

For a short period of time near t ex the bands are in- 
distinguishable in phase space and the oscillation is effec- 
tively period-1. It is this exchange phenomenon that pro- 
duces loops that fill the gap between the period-2 bands 
in the long-time local trajectories (cf. Fig. ||) and con- 
tribute to a sparsely scattered "gas" -like density in G (cf. 
Fig. |). 

An examination of the loop exchanges at different loca- 
tions in the medium revealed the existence of the follow- 
ing spatio-temporal pattern. At any fixed location the 
exchange occurs periodically, with period T ex w 55T r , 
independent of the position (r, 9) in the medium. For 
sufficiently large radii (r > 35) this periodicity takes an 
even stronger form: the entire oscillation pattern, how- 
ever complex, returns with period T ex to the same con- 
figuration. 




X 

FIG. 8. Sketch of the il curve for the disk-shaped array 
(k2 = 1.567, R = 80). Points where the period-2 band ex- 
change was observed are indicated by diamonds. 
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This property smoothly disappears as the defect is 
approached. For two locations rj = (r ,8i) and r 2 = 
(ro, 82) at the same radius ro from the defect but at dif- 
ferent angles, the oscillation pattern at one of them, say 
r2, can be obtained from the corresponding pattern at 
i"i through translation in time by T ex (d2 — #i)/27r, the 
sign of the translation being defined by sign(6>2 — 61). In 
view of this observation it is convenient to introduce a 
coordinate system (r',8') rotating with angular velocity 
2n/T ex relative to the laboratory-fixed coordinate system 
(r, 8). In this rotating frame the local dynamics is de- 
scribed by a time-homogeneous pattern, unique for every 
spatial point r', and the locations in the medium where 
loop exchange occurs correspond to points where the lo- 
cal dynamics always has a period- 1-like character. The 
set of loop exchange points constitute a curve 17 with spi- 
ral symmetry which winds twice around the defect (see 
Fig. H). The two convolutions of 17 lie close to circular 
arcs with radii 19 and 32. This result may be compared 
with the data obtained from an examination of G (cf. 
Fig ^|). The crossings of the bands of G occur at loci 
lying on 17. 

Close to the defect the resolution of the loop exchange 




FIG. 9. Local trajectories calculated on circle with r = 55 for 
scale. 



event is difficult. At r < 18 the difference between the 
period-2 bands is comparable to the band thickness and 
the determination of 17 for smaller radii becomes imprac- 
tical. Variation of the system parameters results in a 
change of the characteristics of 17; for example, the ra- 
dius of the domain R does not affect the shape of the 
17 but does change the angular velocity with which the 
coordinate frame (r',8') in which 17 is immobile rotates 
relative to the laboratory-fixed frame (r, 8) . The angular 
velocity is higher for smaller system sizes: a decrease in 
R from 80 to 60 reduces the period T ex by a factor of 
0.42. A change in the rate constants leads to a defor- 
mation of 17, although the identification of 17 as a set of 
exchange points remains and it retains the topology of a 
curve passing from the defect to the boundary. In Sec. 
we shall show that the existence of 17 is essential for the 
maintenance of spatial continuity in media composed of 
2™-periodic oscillators. 

Simulations on a square array with dimension 80 x 80 
(all parameters were the same as for the disk) show that 
a rotating frame is not necessary to observe the time- 
homogeneous local dynamics of C(r|to,T) = C(y\t). For 
this system geometry the 17 curve is fixed in the medium, 




the square array. All the trajectories are shown on the same 
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a slight wobbling of the defect (frame origin) being ne- 
glected. Figure y shows a number of long-time (r 3> T r ) 
local trajectories on a circle with radius tq = 55 sur- 
rounding the defect in the square domain. One sees a 
significant dependence of the shape of C{tq,Q\t) on the 
angle 9. The local trajectories range from a period-1 orbit 
at the intersection with Vi to the well-established period- 
4 orbit observed in a certain range of 9. To highlight 
the loop exchange phenomenon, a particular time instant 
t = t* is marked on all the trajectories (see Fig. ||) . Com- 
pare the two C(rQ,9\r) at the locations 9\ and 9 2 chosen 
symmetrically on either side of the point 9 = 9^ where 
the circle intersects f2. Visual inspection of these orbits 
shows that their shapes are essentially identical but rep- 
resentative points c(9\,t*) and c((92,i*) appear on dif- 
ferent period-2 bands of the corresponding orbits. This 
clearly demonstrates that the period-2 bands do not just 
approach but indeed pass each other at 9 = 9q , exchang- 
ing their positions in phase space. Since it is not neces- 
sary to work in a rotating coordinate system in the case 
of a square domain, one may resolve the fine structure of 
the local trajectories to a greater degree as can be seen 
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FIG. 10. Cumulative first return map G(0) for the square 
array («2 = 1.567, L = 80, 8 = 0) (top panel) and a magnifi- 
cation of a portion of its structure (bottom panel). Letters on 
the bottom panel denote radii for which corresponding por- 
tions of G(ff) are constructed: (a) 9; (b) 19; (c) 25, and (d) 
31. 



in Fig. [To] (a,b) which shows the cumulative first return 
map G{9) and a magnification of a portion of its structure 
(compare with Fig. ||). The results show that G is com- 
prised of four branches with the fine structure of period-4 
resolved even in the vicinity of the defect (r = 5 is the 
closest distance to the defect for which g(r) is shown). 
Any perturbation of the self-organized pattern of local 
oscillator synchronization due to irregular motion of the 
defect, influence of the boundary or the presence of an- 
other defect may obliterate subtle fine structure of the 
local trajectories. In such a circumstance one is able 
to observe only two gross branches of G and their split 
nature is not resolvable except for very large r. These 
observations allow one to suppose that the local trajec- 
tories may have the same number of fine structure levels 
everywhere in the medium but the degree to which dif- 
ferent levels are resolved in their phase portraits strongly 
depends on the position in the medium relative to the de- 
fect. In view of this hypothesis the phenomenon of spatial 
period-doubling should not be understood in the literal 
sense but rather as an enhanced ability to resolve the fine 
structure with the increase of separation from the defect. 

The stationary rotating spiral wave arises from the 
complex defect-organized cooperation of local oscilla- 
tors. Each location in the medium develops some site- 
specific pattern of oscillation which often differs signifi- 
cantly from that of the corresponding rate law attractor 
and varies substantially from one space point to another. 
There exists a (possibly rotating) reference frame (r', 9'), 
centered on the moving defect, in which local dynamics 
takes a simple, time-homogeneous form. Each point of 
the medium in this frame can be assigned a unique oscil- 
latory pattern, different for different spatial points. This 
allows one to introduce the notion of a defect-organized 
field associated with (r',6') which specifies the pattern 
of dynamics in every spatial point of the medium. This 
field exhibits a complicated architecture lacking of any 
simple symmetries (which can be easily seen from the 
shape of the curve). The slow rotation of this field 
in disk-shaped arrays restores the circular symmetry of 
the solution. Although the manner in which different 
types of local dynamics are distributed in the medium is 
complex, it is not disordered. Due to the continuity of 
the medium maintained by the diffusion, it obeys certain 
topological principles studied in the subsequent sections. 



IV. COARSE-GRAINED DESCRIPTION OF 
LOCAL TRAJECTORIES 



In the previous section the phase space shapes of the 
local trajectories were shown to vary considerably but 
smoothly from one point in the medium to another. To 
describe the transformations of these orbits into each 
other, it is useful to introduce a description which cap- 
tures only topologically significant changes of phase por- 
traits and disregards unimportant details. To understand 
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the topological principles which determine the global or- 
ganization of the defect-organized field one also needs a 
means to compare the time dependence of local trajecto- 
ries. In this section we present a scheme that allows one 
to partition the continuum of all the observed local tra- 
jectories into a finite number of discrete classes according 
to their phase space shape and time dependence. 



A. Representation of attractors by closed braids 

Consider a period-2 rl attractor, P 2 ™, consisting of 2™ 
loops in the concentration phase space V = (c x , c y , c z ). 
Using the cylindrical coordinate system introduced ear- 
lier, we may project P 2 ™ on the (p, 4>) plane preserving 
its original orientation and 3D character by explicitly in- 
dicating whether self-intersections correspond to over or 
under crossings. Such a projection shows a span of </> free 
from crossings where loops are essentially parallel to each 
other. This span can be used to number loops, say, in the 
order of their separation from the origin. This procedure 
maps P 2 " onto a closed braid -E> 2 ™ jig] . Figure [ll] illus- 
trates the construction of the braid representation for the 
P4 attractor of the WR model. 
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It is convenient to subdivide the closed braid P 2 ™ into 
the open braid P 2 »> (separated by dashed lines in Fig. [ll]) 
and its closure where threads run parallel to each other. 
The direction of the flow on the attractor is indicated by 
the arrows. Each crossing on the projection of P4 cor- 
responds to an elementary braid ex.; which refers to the 
fact that thread i overcrosses thread i+1 (cf. Fig. [l7| in 
Appendix for notation rule). An under crossing will be 
designated by o~ x . A braid may be described by a braid 
word that gives the order and types of crossings of braid 
threads. For example, for the closed braid corresponding 
to P4 (cf. Fig. 11) P4 1— > B4 = (T3cr 2 CTif73tT 2 . The closed 
braidP 2 ™ corresponding to P 2 « can be represented by 
several braid words, which can be transformed into one 
another by a set of allowed moves (see Appendix). 

Any braid word representing P 2 n induces a permu- 

(n) 

tation 7T- describing the order in which loops of P 2 n 
are visited during one oscillation period T 2 n . In general, 
each P 2 r. attractor is represented by several possible 7ij , 
their number growing with n; for example, for P 2 there is 

only one permutation tt£ = ( 2 ^) while two permutations 



(2) 



(3421) (which corresponds to the braid shown in 



"1 

Fig. 0) and tt { 2 2) = Qixt) exist for P 4 . With a given 
loop numbering convention each braid word represents 
a unique permutation while one permutation can be in- 
duced by many braid words. 



B. Symbolic representation of periodic orbits 

Take two period-2 n oscillators whose trajectories 
Ci(i),c 2 (i) lie on the same attractor, but which are nev- 
ertheless non-identical since at any given time t their 
dynamical variables are different Ci(t) ^ c 2 (t). Since 
the orbits are periodic there is a time St such that 
Ci(t + St) = c 2 (t) for any t. This operation can be for- 
mally considered as an action of translation operator T$ t 
on the trajectory of the first oscillator: 



T 5t ci(t) = Cl (t + St) = c 2 (t). 



(8) 



FIG. 11. Projection of the P4 attractor on the (c x , c y ) plane 
(top panel) and the corresponding closed braid B4 (bottom 
panel) . 



The concentration time series c(t) of the first oscillator 
then appears to be shifted backward by St relative to that 
of the time series of the second oscillator if St > and 
forward otherwise. Of course, trajectories corresponding 
to different attractors cannot be made to correspond by 
such time translations, e.g. P 2 »> attractors described by 

in] 

different permutations ir\ have different patterns of os- 

(n) 

dilation, but even if two P 2 ™ lie in the same 7r] class 
their actual shapes in V may differ significantly. 

To compare the local dynamics at different points in 
the medium one needs to single out the most important 
characteristic features of the oscillation pattern while 
discarding unnecessary details. A coarse-grained sym- 
bolic description of trajectories appears to be useful for 
this purpose. We assume that the times ti,tn, . . . ,t 2 n at 
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which the trajectory crosses a surface of section = 4>q 
(see Sec. 2) are approximately equally spaced, indepen- 
dent of the choice of 4>q. Thus, the phase point c(t) mov- 
ing along P2™ takes approximately the same time T^n /2™ 
to traverse each loop of the attractor. At t = to 
let the phase point of the period-2™ orbit be on the jo- 
th loop of P 2 n, at t = t + T 2 >./2™ on the ji-th loop, 
and so on (where ji G [1,2™], I G [1,2™], j j^l) until at 
t = to + T 2 ™ the phase point returns to the jo-th loop and 
the pattern (jo, Ji, ■ ■ • J2") repeats. The symbolic string 
s j = (iojii; ■ ■ ■ j2>0 constructed in this way captures the 
most significant gross features of the oscillation pattern it 
describes. In this coarse-grained representation the num- 
ber of possible non-identical trajectories corresponding to 

in) 

a particular 7r 2 of P2™ is finite and the different trajec- 
tories are simply given by the 2™ cyclic permutations of 
Sj. Likewise the time translation operators constitute 
a finite group T h I G [~2 n ~\ 2 n ~ 1 ). They act on the 
symbolic string representing the orbit to give one of its 
cyclic permutations. From its definition it can be easily 

seen that ir^ serves as a symbolic permutation repre- 
sentation of 7+i for the corresponding i-th permutation 
class of Pin . Indeed, consider as an example a period- 
4 oscillation whose representative point lies on loop 3 

at the reference moment of time t — to- Then for the 

(2) 

pattern of oscillation determined by ir\ the state reads 

si = (3241). To obtain the new state translated by T4/4 

backward one acts on si by the permutation representa- 
(2) 

tion ir\ of the 7+i operator to get 

7+1 Si = ( 3241 ) = ( 2413 ) = S 2, 



(9) 



which correctly describes the result of the shift of the 
initial state si. 



V. GLOBAL ORGANIZATION OF MEDIUM 

A. Period-1 regime 

We now return to the spatially distributed medium and 
begin by reviewing some properties of the local dynamics 
in the vicinity of a stable defect with topological charge 
n t = ±1 in a period-1 oscillatory medium. Consider a 
cyclic path r = {r = ?*o > r rnax ,9 G [0,27r)} surround- 
ing the defect. Here r max is a radius such that for all 
(r,6),r > rmax, 9 G [0, 2tt) the shape of the local or- 
bit in phase space V is independent of (r, 9) and closely 
approximates that of the_period-l attractor of the mass 
action rate law (see Sec. |l|). If one starts at an arbitrary 
point (ro,#o) G T one finds that the instantaneous local 
phase 4>(r,t) changes by 2tt or — 2ir (depending on the 
sign of the topological charge) along T. Let us now fix 
a particular time instant t = t* and construct the set of 
points S = {c(r,t*),r G T} as a phase space image of 
instantaneous concentrations at points lying on T. 




FIG. 12. <S-curves (shown by diamonds) constructed for V 
with ro = 55 in period-1 oscillatory (k2 = 1.420) (a) and 
chaotic (k2 = 1.567) (b) media. Solid curves represent short 
time local trajectories on F. 

The property of a defect ([|) and the continuity of the 
medium insure that S is a simple closed curve winding 
once around c*. Figure |l2|(a) shows the 5-curve con- 
structed for the contour T with radius ro = 55, ro > r max 
in a period-1 oscillatory medium with k 2 = 1.420. Since 
all the points on the iS-curve lie at the same time on the 
local trajectories C(r|i ,T), r G T with t* G [to, to + t], 
and for T with ro > r max ah the local trajectories are the 
same and approximated by the period-1 attractor of the 
system (0) , the 5-curve simply coincides with this attrac- 
tor for any t* (cf. Fig. |l2](a)). The iS-curve constructed 
for an arbitrary simple closed path encircling the defect 
in the medium possesses the same property as long as the 
path lies in the open region r > r max . 

This result can be reformulated in terms of time trans- 
lations of local trajectories as follows. Let the local tra- 
jectory C(ro,9o\to,T) at the point (ro,0o) G T be taken 
as a reference, then all of the local trajectories on T can 
be obtained through the translation of C{tq, 6o\to, t) by 
some time St(0 — ) (see Sec. [V). The condition (|) 
implies that 8t{9 — 9q) is a monotonically increasing (de- 
creasing) function such that 6t(2n) = ±Ti where T\ is 
the period of oscillation and the sign is that of n t . Thus, 
the oscillation pattern is continuously time shifted along 
r such that upon return to the initial point it has expe- 
rienced translation by the period. 







B. Period-2™ regime 

For 2™-periodic and chaotic media property (||) holds 
where </>(r, t) should be understood as the angle variable 
introduced in Sec. [y. This can be seen from the following 
argument. Take a period-2 medium with rate constants 
chosen in the vicinity of the bifurcation from period- 1 to 
period-2 such that the attractor P 2 of (J2J) lies infinitesi- 
mally close to P\ from which it bifurcated. Due to the 
continuity of the solutions of the reaction-diffusion equa- 
tion (|3|), the value of § <j)(r, t) dl cannot change abruptly 
when the bifurcation parameter is changed through the 
period-doubling bifurcation. This implies that the S- 
curve constructed for a contour T in a period-2™ medium, 
as in case of a simple period- 1 medium, is a closed curve 
which loops once around c* in phase space. This is illus- 
trated in panel (b) of Fig. [l^ which shows S for contour 
r with radius ro = 55 in medium with = 1.567 and 
time t = t*. Recall again that the points of the iS curve 
have to lie on the local trajectories C(rjto, t), r G T (cf. 
Fig. H where points designated by diamonds lie on S for 
the chosen time moment and contour shown in the fig- 
ure). Since the local trajectories in a period-2™ medium 
loop several times around c*, the curve S which winds 
only once (n t = ±1) around c* cannot span the entire 
local trajectory as is the case for a period-1 medium. As 
one sees from Fig. [l2](b) S follows the larger loop of the 
local trajectory, which for T with ro = 55 is typically a 
period-2 orbit (cf. Fig. ^), and instead of making the 
second turn on the smaller loop, it crosses the gap be- 
tween the loops and closes on itself. Although the shape 
of S changes with time (see for details), for any t* 
there exist segments of S which connect different loops 
of local trajectories. This behavior of the S curves would 
be impossible if loop exchanges were nonexistent. The 
analysis shows that the segments of S covering the gaps 
between the loops of the local trajectories are images of 
points on T which lie close to the intersection with the f2 
curve. Thus, the loop exchanges observed in period-2™ 
media are necessary to reconcile the contradiction be- 
tween the one-loop topology of the S curves determined 
by the presence of a defect and the multi-loop topology 
of the local trajectories determined by the local rate law. 

The change of the local trajectories along the contour 
r in period-2™ media can be considered in terms of time 
translations if one adopts a generalization of the transla- 
tion operation in the following way. In a period-2 medium 
let the contour T and the reference point (r ,6>o) G V 
be chosen so that T intersects the f2 curve in the sin- 
gle point (ro,#n) and suppose that these points are suf- 
ficiently separated from each other. Since the shapes of 
the local orbits change significantly along any closed path 
surrounding a defect (cf. Fig. |^) these trajectories cannot 
be made to coincide by time translation as this operation 
is defined in Sec. IV. Nevertheless, the general features 
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FIG. 13. Period-2 local concentration time series c x (r,t) 
calculated on a cyclic path T surrounding the defect: (a) series 
sampled at four consecutive locations separated by 89 = 30°; 
(b) two series sampled at locations chosen symmetrically on 
either side of the intersection with the Q. curve. 

locations (ro,#i) and (ro,#2) one is able to find a time 
shift At(9i, O2) such that some measure of the deviation 
between the trajectories, say, 



M(At(e 1 ,e 2 )) 



t +T 



z^(t + At)-c^(t)\dt, (10) 



of the temporal pattern of the trajectories are preserved 
(e.g. sharp maxima in Ci(t) time series) and for two 



is minimized. Choosing the local trajectory 

C(r , 8 \to, t) as a reference and comparing it to all 
the other local orbits on Y one is able to define the time 
shift function 5t{9 - 9 ) = At(6,0 o ). The shift function 
St(9 — 9 ) increases (or decreases) monotonically and al- 
most linearly (see Fig. ||(a)) with d{8t)/d9 w T 2 /2 • 2tt 
everywhere on T except for a small neighborhood of 
9 = 9q where it exhibits break. Indeed, the loop ex- 
change at 9 — 9q causes the discontinuity of St(9 — 9q). 
At 9 = 9q both loops of the local orbit become equiva- 
lent and the oscillation is effectively period-1 with period 
T\ = T2/2. Since the loops exchange at 9 = 9q, to find 
the best match ( J10| ) between local trajectories sampled at 
points (r , 9n ~ e) and (ro, #n +e), one needs to translate 
one of the trajectories by St = Ti + 0(e). This can be 
easily seen in Fig. |l^(b) which displays two c x {t) series 
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calculated at spatial points lying 9 - 
side of 9q on T. 



■9n = ±10° on either 



exchanges and allows one to predict the number and the 
kind of loop exchanges necessary to perform the required 
T_ rl . translation. 



C. Trajectory transformations along T 

The transformation of local trajectories along T can 
be imagined to occur a result of two separate processes. 
Suppose everywhere on T except 9 = 9q the shape of the 
local trajectories in V is the same and is equivalent to 
that of C (ro , 9a\to , t) . Then all the other local trajec- 
tories C(ro, 9\to, t), 9 6 [0, 2ir),9 ^ 9q can be found 
by time translation of C(r , 9o\to, r) by St(9 — 9 ) — 
Ti{0 — 9o)/2 ■ 2ir. Assume that all the deformations of 
the phase space portrait of the local trajectory which 
take place along T, including the exchange of loops, oc- 
cur at the point 9 = 9q so that the passage through 6q 
shifts the oscillation by 8t p = T\ — T2/2. Then the result 
of the continuous time translation that occurs during 2n 
circulation along T may be described by the action of the 
% lt operator (n t = ±1), while the result of the loop ex- 
change is described by the operator 71„ t . |l8| The total 
transformation of the local oscillation after a complete 
cycle over T is equivalent to the identity transformation 
and thus the result is in accord with continuity of the 
medium. If one makes the assumption that loop exchange 
does not occur on some contour T encircling a defect with 
\n t \ = 1, the time shift function 5t{9 — 9o) becomes mono- 
tonic and continuous everywhere on T. As a result one 
arrives at the incorrect conclusion that starting from the 
point (ro,9o) with the oscillation pattern symbolically 
represented by the string s±, say s\ — (12), and moving 
along r in the clockwise direction one returns to the same 
point (r , 9q + 2tt) = (r , #0) but with the oscillation pat- 
tern shifted by T2/2 and given by S2 — (21) ^ s\. Note 
that this contradiction does not arise in the period- 1 os- 
cillatory medium where circulation over any closed path 
encircling a defect results in the translation by the entire 
period which automatically satisfies the continuity prin- 
ciple. Thus the necessity of loop exchanges in period-2™, 
n > media with a topological defect demonstrated ear- 
lier in this section in terms of iS-curves is now explained 
in terms of time translations. 

The results for the period-2 medium can be general- 
ized for any n > 1 using the following hypothesis. From 
the main property of a topological defect (|]) it follows 
that integration of an infinitesimal continuous shift d(6t) 
over any closed path surrounding a defect results in a to- 
tal shift by ±T2^ /2" and can be symbolically described 
by the T nt operator. Numerical simulations demonstrate 
the existence of time translation discontinuity points such 
that sum of 5t jumps over these points amounts to a shift 
of ^T 2 ™/2" described by the 71„ f operator. The loca- 
tions of these points in the medium can be identified with 
the f2 curve and the origin of the time translation discon- 
tinuities with the loop exchange phenomenon. The rela- 
tion (|A6|) of the Appendix connects translations and loop 



D. Examples 

Consider again the square 80 x 80 array with rate con- 
stants corresponding to chaotic regime (k 2 = 1.567). As 
period-4 fine structure is the highest level of local organi- 
zation resolved in the medium, it is sufficient to use the 
formalism developed above for P4 to describe the local 
dynamics. The analysis shows that in the bulk of the 
medium the oscillation is given by the tt^ = (3421) pat- 
tern, [fill Using this data and the results presented in 



the Appendix one can easily enumerate all the sequences 
of exchanges resulting in T + \ translation. Indeed, one 
should expect either exchange of loops 3 and 4 followed by 
the exchange of period-2 bands (3412) or first the period- 
2 bands exchange followed by exchange of loops 1 and 2. 
Figure 14 is a schematic representation of the medium 
with a negatively charged (n t = — 1) defect in the center 
and the ft curve displayed. 

Consider the change of the oscillation pattern along 
ray ABC emanating from the defect as the value of 
r increases (see Fig. [w] for the cumulative first-return 
map constructed for this ray) . The pattern of oscillation 

sa = (4132) corresponding to permutation = (3421) 
can be followed from r = 5 to r = 19 where the period-2 
bands undergo exchange. This results in the switch to 

the oscillation pattern described by wjp = Q^ff ) seen at 
(2) 

r = 21. The pattern tt\ ' is restored after loops 1 and 2 
exchange at r — 22 and and this pattern persists until 
another exchange occurs at r = 28. 
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FIG. 14. Sketch of the fl curve for the square array 
(k2 = 1-567, L — 80). The points were obtained from simu- 
lations. The ray ABC intersects Q at locations with radii 20 
and 31. 
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Using translation operator 7+i one can express the 
transition of the state sa {r < 20) through the se- 
quence of loop exchanges described above to the state 
sb = (1324) (for 22 < r < 28) as s B = T +1 s A . The same 
shift can be achieved by continuous translation along the 
path ADEFB which does not intersect fi but winds once 
counter-clockwise around the defect. The c x {t) time se- 
ries at points A, D, E. F and B are displayed in Fig. |l5| 
and demonstrate that this is the case. 




FIG. 15. Concentration time series c x (r,t) calculated at 
points A, D, E, F, B of the square array shown in Fig. 15. 

Continuing to advance along the ray ABC, one finds 
that at r = 28 loops 3 and 4 exchange and oscillation 

(2) 

switches once more to the state corresponding to ir 2 ■ 
After the period-2 band exchange at r = 32 the pat- 

(2) 

tern corresponding to 7rJ is reinstated and remains un- 
changed for all r > 32. Again the oscillation at r > 32, 
described symbolically by sc — (3241), appears to be 
shifted by T4/4 relative to sg and by T4/2 relative to sa- 
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FIG. 16. Segment of the concentration time series c x (r,t) 
calculated for the disk-shaped array (K2 = 1.567, R = 80) at 
r = 76 showing the T4/4 time shift of the oscillation pattern 
(see explanation in the text). 



The existence of a T4/4 shift after crossing f2 can also 
be seen from the results for the disk-shaped array with 
R = 80. Figure |l^ shows a segment of the c x (r, t) time se- 
ries sampled in a fixed frame (r, 9) at r = 76. In this coor- 
dinate system SI slowly rotates clockwise (again n t = —1) 
with period T ex . Two time windows each of length T4 
marked by dotted lines and separated by At = 8T4 allow 
one to see how the oscillation state (4132) is substituted 
by its forward T4/4 translation (2413) after the f2 curve 
passes the observation point at t = t ex . 



VI. CONCLUSIONS 



General principles underlie the organization of 2 n - 
periodic or chaotic media supporting spiral waves. As 
in simple oscillatory media, the core of a spiral is a topo- 
logical defect which acts as an organizing center deter- 
mining dynamics in its vicinity; however, the structural 
organization of the medium that arises from the exis- 
tence of the defect is far more complicated. Due to the 
absence of a conventional definition of phase for oscil- 
lations more complex than period- 1, the identification 
of a defect in terms of the relation (^) is not obvious 
and requires the introduction of (often model-dependent) 
phase substitutes which for some systems may be pro- 
vided by angle variables. Despite the complications with 
the definition of phase, one can identify a defect in terms 
of local trajectories. Indeed, as one moves away from 
the defect the local dynamics takes the form of a pro- 
gression of period-doubled orbits, from near harmonic, 
small-amplitude, period-1 orbits to "noisy" period-2' or- 
bits, where I is a function of variables such as diffusive 
coupling and the system size and shape. 

The presence of a defect imposes topological con- 
straints on the global organization of medium as well. 
As was shown above, when 2" > n t the 2"-loop struc- 
ture of the local trajectories conflicts with the period-n t 
structure of the 5-curve and a complex, asymmetric, spa- 
tial pattern of local dynamics, the defect-organized field, 
arises as a result of the necessity to maintain the conti- 
nuity of the medium. The most prominent characteristic 
feature of this field is the curve defined as the set of 
points where the local dynamics most closely resembles 
period-1. This signals the exchange of period-2 bands. If 
the local trajectories possess structure finer than period- 
2, other loop exchanges leading to more subtle changes 
in the local orbits can be found in the vicinity of f2. The 
net result of these exchanges is to produce a time shift of 
the trajectories which compensates for the smooth time 
translation accumulated on continuous paths. Since the 
topological continuity must be observed on any arbitrar- 
ily large closed path encircling a defect and, therefore, 
this contour has a point of intersection with f2, a single 
defect in a period-2™ (n > 0) medium cannot be local- 
ized. 
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We point out again that many of the phenomena we 
have discussed above are not dependent on the existence 
of a period-doubling cascade or chaotic local dynamics, 
although this is the case we have analysed in detail. 
Reaction-diffusion systems with local complex periodic 
orbits in phase space dimensions higher than two should 
exhibit similar features when they support spiral waves. 
It should be possible to experimentally probe the phe- 
nomena described in this paper. The appropriate pa- 
rameter regime can be determined from investigations of 
well-stirred systems. For example, period-doubling and 
chaotic attractors have been observed in the Belousov- 
Zhabotinsky reaction. [ p0[ If the spiral wave dynamics is 
then studied in a continuously-fed-unstirred reactor [^l) , 
one should be able to observe the characteristics of the 
spiral dynamics and the loop exchange process that serve 
as signatures of the phenomena described above. 
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APPENDIX A: BRAID MOVES AND LOOP 
EXCHANGE OPERATORS 

In this appendix we make use of the projection of the 
period-doubled attractors P2™ onto closed braids £?2« (see 
Sec. IV) to show how loop exchanges affect the pattern of 
We demonstrate that those combinations of 



oscillation. 

loop exchanges that produce identity transformations of 
Pin result in nontrivial time translations of trajectories. 

Each closed braid is represented by a set of non- 
identical braid words with their number rapidly growing 
with n. Without violation of the topology of they 
can be transformed into one another by the following set 
of moves (see, e.g. ]l6|) : 

1. commutation relation, <Ji<7j = <Jj<Ji, \i — j\ > 2; 

2. type 2 Reidemeister move, <Ji<J^ — cr^~ ai = 1; 

3. type 3 Reidemeister move, OiOi+\Oi — Oi+iOi<Ti+i, 



4. first Markov move, <7;E<7, ; 



E, E G B; 



where B is a set of open braids. While the first three rules 
are common for all braids, rule 4 is specific for closed 
braids. Indeed, it can be written in the form <7;E = Ecr.; 
which, for elementary braids o"j, corresponds to moving 
Cj around E>2 n resulting in the exchange of the closed 
braid loops (cf. Fig. |l7|(d)). Type 1 Reidemeister (or 
second Markov) moves are not allowed since they do not 
preserve the number of loops, an essential feature of P2™ 



FIG. 17. Conventional designations and basic braid moves: 
(a) definition of elementary braids Oi and er" 1 ; (b) type 2 Rei- 
demeister move; (c) type 3 Reidemeister move, and (d) first 
Markov move (E represents arbitrary braid). 

attractors. While rules 1, 2 and 3 do not affect the 
first Markov move does (except for the degenerate case of 
P2 which is represented by the single permutation tt^). 
Thus any number of rearrangements affecting only the 
braid -B2™ of the closed braid £?2™ leave the braid word 
in the same permutation class 7r| , while each applica- 
tion of the first Markov move yields a new permutation 
class. 



1. Loop exchanges for P2 and P4 



We now examine how the loop exchanges influence 
the patterns of oscillation for the period-2 and period- 
4 attractors. For Pi one has only the single braid word 
a'i((Tj _1 ) and the single permutation 7^ = (21) induced 
by (7 1 . Two different symbolic states s\ = (12) and 
S2 = (21) are possible for the period-2 oscillation with 
respect to some fixed time frame. We introduce an oper- 
ator whose action on the closed braid representing 
P2 is to move a\ by 2-7T in a direction opposite to the 
flow. The result of the action of this operator, which is 
the first Markov move for E = 1, is to leave the attractor 
P2 unchanged; however, one finds that loops 1 and 2 have 
exchanged their locations in phase space. In the time se- 
ries for the dynamical variable Ci (t) the exchange can be 
seen as a substitution of taller maxima (2) by shorter 
maxima (1) and vice- versa. If this process is followed 
in time it produces the characteristic pattern shown in 
Fig. 0. Thus, application of A± to P 2 induces a transfor- 
mation of the oscillation state s\ into S2 and vice- versa. 
This can be symbolically described as an action of an 
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exchange operator A^ represented by the permutation 
(21) : 



4 1] *i = (21) (12) = (21) = s 2 , 



4^ *2 



(if) (* 



(121 



One sees that the action of 4 is equivalent to that 
of 7+i which translates the oscillation pattern by half a 
period. The inverse of the braid operator A± can be in- 
troduced in an analogous way as an operator moving a\ 
along the direction of the flow. It corresponds to an ex- 
change operator (4 ) _1 = 71 1 acting on strings. Since 
for P 2 application of T +1 or T_i results in essentially the 
same states the sign of the shift is chosen to maintain 
consistency with corresponding operators for P 2 n, n> 1. 

Double application of 4^ results in translation by a full 
period and thus in the identity operator 



Si. 



(Al) 



(A\ 



(lh2 



(A\ 



0-h 



(A2) 



The P4 attractor possesses a richer set of transforma- 
tions. Since under coarse-graining braid words which in- 
duce the same are indistinguishable, we can single 

(2) 

out two essential representatives Si = 03 S for n\ — 

Qfjf) and S 2 = (7i S for ir {2) = Q|f|), where S stands 

(2) - 
for <72cri(73f72- Let A\ be an operator on B 2 which 

moves the double-thread crossing (cf. large dashed box 

in Fig. |ll|) S by 2n in the direction opposite to the 

flow, analogous to the action of 4 on o\ for B 2 . In 

(2) 

fact, the action of A\ ' can be seen as an exchange of 
the period-2 bands of P4, each consisting of two period- 
4 loops. Unlike 4 the operator 4 alternates braid 
words and the corresponding pattern-defining permuta- 



tions £i,7rj 2 ' ) <— — > ^,11^'. The application of Af in- 
duces transformations of symbolic strings Sj which again 
can be described by the action of an exchange operator 
1(2) 



J2) 



1 (2) 



the string transformation described by the exchange op- 



erator A 2 2 ^ with symbolic representation (J^lf)- When it 
acts on E2 = (7i S it exchanges loops 1 and 2 and the per- 
mutation representation of 4^ changes to (2134)- The 

(2) 

inverse of A 2 ' moves the single crossing along the flow 
and produces opposite results; i.e., acting on Si it leads 
to 1 o 2 exchange and when applied to S2 it results in 
the exchange 3^4. 

Note the difference between action of 4 on braids 
and the action of exchange operators 4™'* on symbolic 
strings. While several operations 4 applied to the 
same initial braid word lead to equivalent final words 



(e.g. A^ and A 2 Z> ) the resulting loop exchanges and, 
thus their permutation descriptions, can be quite differ- 
ent ((3412) and (1243) for the example chosen). Conse- 
quently, compositions of braid operations returning the 



,(2) 



braid B 2 n to its initial state (e.g. 
duce nontrivial translations of Sj. 

let {A 2 2 ^) 2 act on the trial state si 



(2) (2) 
A 2 ' o A 2 ) may in- 

To demonstrate this 
= (3241). Applying 



the rules one obtains 

(2134) (l243j (3241) = (4132) = r+ 2 .si, 

(A3) 



(4 2) ) 2 si - fi 2 ^ ^ 123 ^ 



thus relating the simultaneous loop exchange (13) <-> (24) 

with translation of the period-4 oscillation by half of a 

(2) 

period. This implies as well that application of A 2 
four times results in the identity string transformation 
(4 2) ) 4 =7+4 = 1 and, therefore, = {A {2) f. 



(2) 

and A 2 provide 



(2) 

Compositions of braid operators A\ ' 
another example of how identity braid operators induce 
nontrivial string transformations. Since both operators 
and their inverses alternate braid words Si <-> S2 the ap- 
plication of the composition of any two of them returns 
(2) 

the braid to the same tt] permutation class and, thus, 
the resulting string transformation is equivalent to some 

translation. The relations for the compositions of the 

(2) (2) 

A\ and A 2 operators can be obtained directly from 



A\ represented by the permutation Due to the their symbo i ic representations 



apparent similarity of action of A± and A\ z> on the cor- 

(2) 

responding attractors, A\ inherits the algebraic prop- 
erties of Ai'. Indeed, 4 produces identity operator 
when applied twice and, thus, is equal to its inverse. 

Finer rearrangement of the P4 loop structure is pro- 

(2) 

vided by the action of A 2 defined as an operator which 
moves the single crossing (enclosed in the smaller box in 
Fig. O) by 2-7T in the direction opposite to flow. From 

the structure of B4 one sees that after application of A 2 ' 
the single crossing does not return to the same location 
in V but appears on the other period-2 band; thus, braid 



,(2) 



A2) A (2) _ H234\ (Vm\ _ A (2) A (2) 

Ai oA 2 - ^243 ) ^3412 ) - A 2 A 



A\ 



(2) 



/1234\ /T234\ 
\M12J (,2134 J 

l(2h-l 



/1234\ 
(,3421 ) 

H234\ /1234\ 



,( 2 ) 



71 



(A4) 



°wr l = [%ffi) (3li2j = (4 2) r 1 °4 2) 



_ H234\ /T234\ _ /T234\ _ fa) 
- \M12J (,1243 J - \A312J ~ n 2 



= 71 



words and n. 



(n) 



permutations alternate. Depending on 
2^ results in differ- 
ent loop exchanges. In the case of Si = 03S action of 

(2) 

A 2 leads to the exchange of loops 3 and 4 and results in 



the initial state, the application of A 9 



These relations are constructed using the assumption 
that the initial state of the braid is Si. Although appli- 
cation to an alternative initial condition changes actual 
permutation representations of the exchange operators it 
yields algebraically equivalent results. From (A.4) one 
sees that all the exchange operators commute and their 
compositions provide operators which translate the oscil- 
lation by all the allowed multiples of T 4 /4. 
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2. Loop exchanges for P2™ attractors 



A generalization of the phenomena discussed above 
to arbitrary n may be inferred from the observation of 
the structural organization of the closed braids i?2™ cor- 
responding to period-doubled attractors P^n. Indeed, 
B 2 n+i can be obtained from B^ by doubling each thread 
of -B2™ and adding a single crossing on top to preserve 
simple connectivity of the construction. The braid i? 2 ™ 
arising as a result of n successive iterations of this pro- 
cedure can be subdivided into n non-overlapping struc- 
turally similar blocks of braids Em', m—l,n. This prin- 
ciple of structural organization is illustrated in Fig. [l8| 
representing Bg and its three crossing blocks shown in a 
series of boxes with decreasing sizes. The analysis shows 
that these blocks can be moved as whole entities along 
i?2" without interference from each other resulting in the 
exchange of those loops along which they move. The es- 
sential parts of these moves can be represented by a set 
if of 2ir movements of structural blocks Em' so that 
corresponds to the largest block and results in an 

exchange involving all the 2™ loops, A 2 corresponds to 
movement of next-smaller braid block and results in the 
exchange of 2 n_1 loops, and so on. The transformations 
of time trajectories resulting from exchanges of loops can 
be again described by the action of permutation opera- 
tors Am on symbolic strings Sj . The fact that the cross- 
ing blocks move independently results in the commutivity 
of operators A^m with each other. The geometry of £?2™ 



also defines the basic algebraic property of A- 
strated above for the n = 1, 2 examples 



(n) 



= 1, m £ [1, n]. 



demon- 



(A5) 



Some compositions of the exchange operators yield trans- 
lation operators % where I £ [— 2 n_1 , 2™~ 1 ). For the dis- 



cussion of the phenomena described in Sec. Ill only the 
operator 7+i and its inverse are of particular interest. 




Using induction from the analysis of cases with small 
n one may infer the general expression for the T + i trans- 
lation operator : 



(A6) 



m— 1 



FIG. 18. Braid B» constructed for the Pg attractor. 
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